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Abstract: Analysis of data from the Rossi X-ray Timing Explorer satellite's All Sky Monitor instru- 
ment for several X-ray binary sources has identified a recurrent ~24 h period. This period is sometimes 
highly significant, giving rise to the possibility of it being identified as an orbital or super-orbital pe- 
riod. Further analysis has revealed the same period in a number of other X-ray sources. As a result 
this period has been discounted as spurious, described variously as arising from daily variations in 
background levels, the scheduling of ASM observations and beating between the sampling period and 
long-term secular trends in the light curves. We present here an analysis of the spurious periods and 
show that the dominant mechanism is in fact spectral leakage of low-frequency power present in the 
light curves. 

Keywords: instrumentation: miscellaneous — methods: data analysis — space vehicles: instruments 
— techniques: miscellaneous — telescopes — X-rays: binaries 
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1 Introduction 

The Rossi X-ray Timing Explorer (RXTE) spacecraft 
was launched on 1995 December 30 into a circular low- 
inclination (~23°), Low Earth Orbit (~580 km), with 
an orbital period of ~96 min (Jahoda et al. 1996). It 
continues to operate almost flawlessly up to the present 
time. The payload consists of two pointed instruments 
- the Proportional Counter Array (PCA) and the High 
Energy X-ray Timing Experiment (HEXTE); and one 
scanning instrument - the All Sky Monitor (ASM). 

Analysis of ASM data using a variety of techniques 
has identified a number of spurious periods. Power at 
the ~96 min spacecraft orbital period is evident in 
the power spectra of most sources, with other spuri- 
ous periods reported at ~180 d and ~365 d (Benlloch 
2004). Additional periods at ~50 d and ~11 yr might 
be expected due to the precession of the spacecraft or- 
bit (Corbet 2003), and from modulations due to the 
solar cycle, although the ASM data sets are not long 
enough to see an ~11 yr periodic variation as any- 
thing other than a long-term secular trend. A highly 
significant period at ~24 h has also been reported for 
a number of sources and has been variously attributed 
to daily variations in the X-ray background (Corbet 
et al. 1999), beating between the sampling period of 
the satellite and a long-term secular trend in the light 
curve (Benlloch 2004), and to the scheduling cycle of 
the ASM (Zdziarski et al. 2004). At least two X-ray 
binary systems are suspected to have orbital periods 
around this value: GX 339-4 (~0.5 - 1.7 d; Buxton & 
Vennes 2003; Hynes et al. 2003) and CAL83 (1.0436 
d; Smale et al. 1988). While variations at the orbital 
period in X-rays do not appear to be evident in RXTE 
ASM data for either source, it is entirely possible that 
other sources exist which exhibit real variations around 
this period. It is thus important that the underlying 



mechanism be determined so that spurious periods can 
be confidently ruled out. 



2 The All Sky Monitor 

The ASM consists of three wide-angle Scanning Shadow 
Cameras (SSCs) covering a nominal 1.5 - 12.0 keV en- 
ergy range with three energy channels (1.5 - 3.0, 3.0 
- 5.0, and 5.0 - 12.0 keV; A. M. Levine 2004, private 
communication). The ASM assembly scans ~80% of 
the sky every orbit, dwelling on a target source for ~90 
s each observation (Bradt, Rothschild & Swank 1993). 
Each camera consists of a position sensitive xenon pro- 
portional counter covered by an 8 /im thick aluminized 
plastic thermal shield and a 50 /im thick beryllium win- 
dow. In addition, the proportional counters in SSCs 2 
and 3 have a ~2 fim thick coating of polyimide on the 
interior of the beryllium foil as a protection against 
leaks (Levine et al. 1996). The field of view of 6° x 
90° FWHM is limited by a thin aluminium coded mask 
collimator. The current target catalogue contains over 
300 X-ray sources, each observed up to 15 times per 
day by one or more of the SSCs. As a result, the ASM 
light curves can contain a number of data points within 
any given orbit, as the source moves in and out of the 
fields of view of the three SSCs. 

The ASM collects data with a duty cycle of ~40% 
which, when combined with spacecraft manoeuvres, 
results in a highly stochastic pattern of sky coverage 
(Levine et al. 1996). The resulting observing pat- 
tern for an individual source is thus uneven, with the 
orbital location of RXTE while observing a specific 
source seeming to change randomly over time. The 
pattern of the shadow cameras coded mask along with 
a model of the diffuse X-ray background is taken into 
account, with the source intensities obtained from the 
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solution of a linear least-squares fit. Corrected inten- 
sities from individual dwells with high values of re- 
duced x 2 , short exposure times due to entry into high- 
background regions, from sources close to or behind 
the Earth's limb, from sources close to the edge of 
the field of view, and from fields that are exception- 
ally crowded are rejected (Levine et al. 1996). The 
resulting background subtraction is not perfect, with 
numerous negative count rates appearing in each light 
curve. The corrected intensities from individual ob- 
servations (or dwells) are compiled and made available 
to the public domain on NASA's High Energy Astro- 
physics Science Archive Research Centre (HEASARC) 
website 1 . The ASM has been operating on a nearly 
continuous schedule, turning off only during passage 
through high-background regions such as the South 
Atlantic Anomaly (SAA) and whenever the Sun is in 
the field of view (Levine et al. 1996). 

Other systematic errors that may affect some sources 
include effects upon the intensity when another bright 
source is in the field of view, incompleteness of the 
active catalogue of X-ray sources, and low energy ab- 
sorption by the window and thermal shield increasing 
with elevation of a source in the field of view. In addi- 
tion, a slow drift in the gain in SSCs 1 and 3 has been 
observed over time, along with a change in the chan- 
nel definitions used to define the three energy bands 
which occurred at MJD 51548 (Corbet 2003). SSC 3 
has sustained the most degradation over time, while 
the gain in SSC 2 has remained relatively stable. 

3 The -24 h Period 

3.1 Data Analysis and Reduction 

RXTE ASM background-subtracted single dwell data 
were used for all our analyses. As the observed 24 h 
periods are not believed to originate from the distant 
source, a barycentric correction was not performed. 
Data from SSC 3 were rejected due to the degrada- 
tion of the instrument, and linear trends were fitted to 
the light curves before and after the gain adjustment 
in MJD 51548. These linear trends were then sub- 
tracted to remove the resulting spurious low-frequency 
power. All individual dwells with an error greater than 
Ict from the mean error in the light curve were re- 
moved, along with any dwells with anomalously low 
count rates. The mean count rate was then subtracted 
from the light curve to reduce noise in the power spec- 
tra. 

The primary analysis technique employed on the 
reduced light curves was the Lomb-Scargle periodogram 
(Lomb 1975; Scargle 1982), although epoch folding and 
epoch folding search techniques using \ 2 statistics were 
also employed to test that the problem was not a re- 
sult of the analysis technique alone. The Lomb-Scargle 
power spectrum is particularly useful in the identifica- 
tion of periodic signals buried in noisy unevenly spaced 
time series, while the epoch folding search technique 
is best suited for the identification of non-sinusoidal 
signals. Epoch folding of a light curve at a suspect 

1 http://heas arc. gsfc.nasa.gov/docs/xte/xtegof.html 



period can then determine the shape and parameters 
of a periodic signal. All three of these techniques de- 
tected the presence of ^24 h modulation in the sus- 
pect light curves, indicating that the spurious period 
is not an artefact of a particular analysis technique. 
For the periodogram analyses we employed the fasper 
subroutine (Press & Rybicki 1989), while the epoch 
folding and epoch folding search techniques made use 
of the FTOOLS analysis software available through the 
HEASARC website. It should be noted that the False 
Alarm Probability (FAP) of a given peak in a Lomb- 
Scargle power spectrum depends on the sampling pat- 
tern of the data (Home & Baliunus 1986). Therefore, 
we used Monte Carlo simulations (e.g. Kong, Charles 
& Kuulkers 1998) to determine the 99% white noise 
significance levels for all power spectral analyses. 

3.2 Characterisation of the ~24 h 
Period 

During the course of analysis of ASM light curves, 
we first identified a spurious period around 24 h for 
the Low Mass X-ray Binary black hole candidate GX 
339-4. This happened while searching for its previ- 
ously reported orbital period of between ~0.5 - 1.7 
d. Highly significant power was observed in the Lomb- 
Scargle power spectra around 24 h, which we provision- 
ally associated with the orbital period. However, fur- 
ther analysis of ASM data revealed four other sources 
within 5° of GX 339-4 that also showed power at ~24 
h. Benlloch (2004) had also noted the 24 h period 
in the GX 339-4 light curve, attributing it to beating 
between the sampling period of the spacecraft and a 
secular trend in the light curve. Zdziarski et al. (2004) 
on the other hand pointed to the scheduling cycle of 
the RXTE observations as the cause of the presence of 
this period. 

Furthermore, Corbet, Finley & Peele (1999) had 
identified a similar 25.65 h period in the High Mass 
X-ray binary pulsar 2S 0114+650, and had explained 
it as being related to daily variations in background 
levels. They had also noted that variability at similar 
periods had been seen in at least some other sources 
observed with the ASM. The observed variability in 
2S 0114+650 was attributed to passage of the space- 
craft through the SAA (R. H. D. Corbet 2004, private 
communication). While this 25.65 h period is clearly 
present in the ~2 yr of ASM data that Corbet, Finley 
& Peele utilised, it is not significantly present in our 
8.5 yr data set (see Farrell, ONeill & Sood 2005). Sub- 
sequently, many other sources have been revealed that 
show power at ~24 h. 

Perhaps the most puzzling aspect of the 24 h period 
is that it is not present in data from all sources - it 
seems to appear only in conjunction with significant 
low-frequency power, although there are cases where 
sources with confirmed low-frequency modulation do 
not show significant power at or around the suspect 
period. The sources with observed 24 h periods are 
not all grouped together in the sky but tend to vary 
in their celestial coordinates. Furthermore, the value 
of the period is not constant between sources, and is 
sometimes seen to split forming a double peak either 
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side of 24 h. The power of the period varies between 
sources and even within a single light curve, sometimes 
disappearing completely in certain sections of a data 
set. 



3.3 Sources with ~24 h Periods 

We have identified 22 X-ray sources with power at ~24 
h in their power spectra. Figure shows the location 
in the sky of these sources, grouped into sources with 
strong, medium, and weak power at 24 h. From these 
we selected six well studied sources randomly located 
in the sky for our analyses, each marked on Figure Q 
and referred to in all subsequent figures with the cor- 
responding letters: Crab pulsar (A), LMC X-4 (B), 
Vel X-l (C), Her X-l (D), GX 339-4 (E), and EXO 
2030+375 (F) . Each of these selected sources showed 
evidence of the 24 h period with > 99% significance in 
all available ASM data (~8 yr of data), and in a 400 
d time period spanning MJD 52385 - 52785. Figure H 
shows the power spectra of all available ASM data for 
each of these sources, covering 0.1 - 1 000 d. While the 
low-frequency power present in the Crab power spec- 
trum is unexpected, the power spectra of the other five 
sources matches well with previously published data, 
indicating that our analysis techniques are sound. The 
peaks present at periods > 150 d in the Crab power 
spectrum are indicative of long-term variations in the 
light curve, possibly due to background variations aris- 
ing from the low declination of the source and/or the 
solar cycle. Figure is a higher resolution look at the 
24 h peaks in the same power spectra covering 0.8 - 
1.2 d, with the 99% significance levels indicated. 
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Figure 1: Location of 22 sources identified as hav- 
ing ~24 h periods, ranked in terms of strong (dia- 
monds), medium (circles) and weak (crosses) sig- 
nal/noise ratio of the 24 h peaks. Sources used for 
our analyses were: (A) Crab pulsar; (B) LMC X-4; 
(C) Vel X-l; (D) Her X-l; (E) GX 339-4; and (F) 
EXO 2030+375 
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Figure 2: Power spectra of the six test sources us- 
ing all available RXTE ASM data. Each source is 
referenced using the same letter denoting its loca- 
tion in Figure n 
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Figure 3: Power spectra of the six test sources us- 
ing all available RXTE ASM data showing the 24 
h peaks. The 99% significance levels are indicated 
by the dashed lines 
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4 Daily Variations in 
Background Levels 

4.1 SAA Passage 

The area of low magnetic rigidity known as the South 
Atlantic Anomaly has long been known to cause se- 
rious problems for spacecraft. The SAA is believed 
to be a result of the magnetic dipole axis being offset 
from the centre of the Earth, resulting in a compara- 
tive weakening of the geo-magnetic field off the coast 
of Brazil. An increase in the radiation levels at Low 
Earth Orbit altitudes can lead to single event upsets, 
surface charging, and component failures among other 
problems. As a result, most spacecraft routinely switch 
off critical subsystems during passage over the SAA in 
order to minimize disruption and damage. 

The high energy of the charged particles in this 
region is often within the range of interest of satel- 
lite based X-ray and gamma-ray telescope instruments. 
The very high fluxes of these particles (> 1 000 pro- 
tons cm~ 2 s _1 above 30 MeV) result in the activation 
of isotopes within the spacecraft material during pas- 
sage (Shaw et al. 2003). These isotopes have vary- 
ing decay half-lives, sometimes lasting until well after 
the spacecraft has left the SAA region. During this 
decay time-scale the background radiation levels are 
increased, resulting in a higher count rate from in- 
struments such as the proportional counters used in 
the ASM and the PCA. Observations made within the 
time frame of activation decay inherently introduce pe- 
riodicities into the light curves such as the spacecraft 
orbital period and the rotational period of the Earth. 
RXTE for instance regularly passes through the SAA 
approximately 10 out of every 15 orbits, as the Earth 
rotates slowly under the orbital path (E. Smith 2004, 
private communication) . Background levels due to ac- 
tivation will gradually increase with each orbit until 
the SAA moves out from underneath the orbital path. 
The background levels then slowly decrease for the 5 
orbits until re-entry, hence introducing an approximate 
24 h periodicity into the data. Another periodicity at 
~50 d may be introduced in a similar fashion as the 
orbit of RXTE gradually precesses at this period. Af- 
ter 50 d the time-of-day that SAA passages occur will 
have moved earlier through one whole day, resulting 
in 51 cycles of SAA passage within the 50 d (E. Smith 
2004, private communication). 

The differing decay half-lives of radioactive iso- 
topes contribute separate components to the background 
levels of both the PCA and ASM instruments. The two 
most significant components in the PCA background 
have ~25 and ~240 min decay time-scales, although 
the specific contributing species have not been iden- 
tified (C. Markwardt 2004, private communication). 
The instruments that make up the PCA consist of five 
xenon/methane proportional counters with a 90/10 
mixture ratio and a propane veto level (Jahoda et al. 
1996). Each ASM SSC consists of a xenon/C0 2 pro- 
portional counter with a 95/5 mixture ratio (Levine 
et al. 1996). As both instruments are xenon propor- 
tional counters, it is fair to assume they will experience 
similar activation decay time-scales. 



The main material components for both instru- 
ments are xenon and aluminium. The background in- 
troduced into the detectors as a result of the passage 
through the SAA will be a function of the relevant 
excitation cross-sections for these elements, the par- 
ticle flux in the SAA, and the half-lives of the radio- 
isotopes consequently produced. Aluminium does not 
have radio-isotopes that have characteristic activation- 
decay time-scales of ~25 or ~240 min. 120Xe, 121Xe 
and 123Xe have half-lives of 40, 39, and 120 min re- 
spectively, but it is highly improbable that these iso- 
topes will be produced in the expected scheme of ex- 
citation, as the most abundant naturally occuring iso- 
topes of xenon have atomic mass numbers of 129 or 
greater. 
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Figure 4: Power Spectra of the 400 d light curves 
spanning MJD 52385 - 52785 for the six test 
sources 



In order to investigate SAA passage as the root 
mechanism for the observed 24 h periods, we first mod- 
elled the location of the SAA with respect to the space- 
craft over time. To do this we made use of the RXTE 
orbital ephemeris data file available via the HEASARC 
website. Due to the high time resolution of these data 
(the orbit file contains the location and velocity of the 
spacecraft every 60 s since launch) and the size of the 
data file, we selected a 400 d section spanning MJD 
52385 - 52785 for our analysis. Using the RXTE house- 
keeping data from PCA observations of 2S 0114+650, 
we were able to determine the location of the SAA 
relative to the spacecraft at one particular instance in 
time. The housekeeping files give the spacecraft loca- 
tion in Earth latitude and longitude against mission 
time, allowing us to calculate the RA and dec of the 
SAA for a given time stamp. Taking the coordinates 
of the SAA maximum as 25° S, 45° W (Plett 1978) 
we were able to model its movement by keeping the 
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dec constant and stepping the RA parameter with the 
23.9345 h rotational period of the Earth. In order to 
validate our model we referred to PCA observations 
of the Crab pulsar, Cir X-l and 2S 0114+650 listing 
the time since SAA passage for each individual time 
increment. We compared the location of the SAA ob- 
tained from our model with that of the spacecraft dur- 
ing small values of the time since SAA parameter, and 
found our model to be a good fit for all three sources 
over the 8.5 yr mission life of RXTE. 
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Figure 5: Power spectra of the 400 d light curves 
after SAA affected observations were filtered out 

We generated power spectra for the six test sources 
using the Lomb-Scargle periodogram technique for the 
reduced ASM data from SSCs 1 & 2, between MJD 
52385 - 52785 (Figure 0}. Applying our model of the 
SAA coordinates, we generated a list of times during 
which RXTE was within a conservative region of ± 
45° RA and ± 25° dec centred on the location of the 
SAA maximum. A new time since SAA passage pa- 
rameter was then calculated for each individual dwell 
time stamp for each of the six test sources. These 
parameters were used to filter out observations that 
were made within a given time (AT) of the spacecraft 
leaving the SAA affected region. Filtering out all ob- 
servations within an initial value of 30 min for AT 
and recalculating the power spectra, we found that 
the 24 h period was not significantly reduced for any 
of the sources. Repeated analysis with a AT of 120 min 
should have seen significant reduction in the power of 
the 24 h peak, due to the removal of all observations 
within one orbit of SAA passage. Filtering with a AT 
of 120 min would also remove all observations possibly 
affected by the decay of any xenon isotopes activated 
by SAA passage. Our results show that this is not the 
case - not only were the peaks still present they had 
significantly increased in power. Figures [I] & show 



the power spectra for each of the test sources before 
and after SAA filtering. While it is possible that the 
decay of an isotope with a longer half-life (such as the 
240 min component mentioned previously) could be 
causing the 24 h period, this is unlikely due to an ad- 
ditional factor. Simulated power spectra of reported 
super-orbital periods present in a number of the af- 
fected sources (e.g. Her X-l & GX 339-4) show power 
around 24 h, indicating that the spurious period is 
most likely not an artefact of the instrument count 
rates but instead related to the sampling period. 

4.2 An ' Anti- Anomaly' ? 

While the offset of the geo-magnetic dipole results in 
an area of low magnetic rigidity at the SAA, an area 
of high magnetic rigidity must exist on the oppos- 
ing hemisphere, which we will refer to as the Anti- 
Anomaly (AA). Passage through the A A would result 
in lower levels of background radiation due to increased 
shielding from the geo-magnetic field. This could re- 
sult in periodic dips in the background levels, possi- 
bly manifesting in the power spectra as peaks simi- 
lar to those described previously for the SAA passage. 
Spacecraft instruments and subsystems are not com- 
monly shut down during AA passage so a number of 
observations are bound to occur over the AA each or- 
bit. The issue of isotope activation is not present in 
this case due to the reduction in the flux of high energy 
charged particles. Taking the location of the AA as 10° 
N, 100° E (Plett 1978), we modelled the movement of 
the AA in a similar fashion to the method previously 
described for the SAA. Using the same six test sources, 
we filtered out all observations made within a region of 
± 45° RA and ± 25° dec centred on the AA location. 
Figure |S| shows the power spectra for each of the test 
sources after AA filtering. The power spectra shown 
in Figure [5] are practically indistinguishable from the 
unfiltered power spectra shown in Figure 0] indicat- 
ing that AA filtering did not significantly reduce the 
power of the 24 h peaks and thus is not the dominant 
mechanism. 

5 Spectral Leakage 
5.1 ASM Sampling Rate 

The ASM observing schedule is determined by a soft- 
ware package that takes into account a number of vari- 
ables, generating observing plans for one day time in- 
tervals. These variables include the pointing program 
of the spacecraft, the spacecraft orbital ephemeris, the 
location of the Sun and the apparent direction of the 
Earth, the limits of the ASM drive assembly rotation 
range, the slew rate of the drive assembly, and the loca- 
tions of high background regions. The daily observing 
plans are highly complex, such that the times of ob- 
servation of a given source tend to look random (A. 
M. Levine 2004, private communication). The sam- 
pling rate for a given source is partially uneven, with 
blocks of unevenly spaced observations reoccurring on 
a quasi-periodic basis. It is possible and even probable 
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Figure 6: Power Spectra of the 400 d light curves 
after AA affected observations were filtered out 



that this grouping has a period of around 24 h. 

In order to test whether the sampling rate of ASM 
observations is the cause of the spurious periods, we 
simulated light curves with the same sampling rate of 
the affected sources but with random intensity values. 
The resultant Lomb-Scargle power spectra for each of 
the sources showed no significant power at any period, 
indicating that the sampling rate alone cannot be the 
cause of spurious periods at 24 h or any other period. 
When the random intensity values were replaced with a 
30 d period sine wave, power at 30 d and ~1 d emerged 
in all the power spectra. Closer examination of the ~1 
d peaks resolved them into two closely spaced peaks. 
Repeating the simulations with different input periods 
found that the spacing between the I d peaks was in- 
versely proportional to the value of the input period, so 
that as the input period was reduced the peaks moved 
apart, and as it was increased the peaks converged. 
These results strongly imply that the spurious 24 h 
periods are not related to variations in background lev- 
els, but instead result from an interaction between the 
sampling rate and real modulations present in the light 
curves. 

5.2 Spectral Leakage Effects 

Spectral leakage is the transferral of power from a cer- 
tain frequency to either higher or lower frequencies in a 
time series. Leakage to nearby frequencies is due to the 
finite length of the time series, while leakage to distant 
frequencies is a result of the finite size of the interval 
between samples (Scargle 1982). Unevenly sampling a 
data set tends to reduce and often eliminate the effect 
of spectral leakage. While astronomical data by the 
nature of the observations are typically unevenly sam- 



pled, there is sometimes an enforced regularity which 
results in semi-regular spacing of data points (Scar- 
gle 1982). Events that can enforce regularity in the 
case of the RXTE ASM include Earth occultation of 
the source, regular shut-down during passage through 
high background regions (e.g. the SAA), and the orbit 
of the Earth around the Sun. This enforced regularity 
can quite often lead to spectral leakage of a signal to 
other frequencies, and in fact explains why filtering of 
all SAA affected dwells increased the power of the 24 h 
peaks. Removal of all observations within 120 min of 
an SAA passage inherently increases the enforced ~24 
h regularity, increasing the leakage of power. 

Benlloch (2004) found that leakage from low to 
high frequencies does occur in ASM light curves owing 
to the sampling of the data. While noting the general 
cause of the 24 h period in GX 339-4, the pseudowin- 
dow function of the data was not examined and the 
effect was not explored in-depth by this author. As 
such the appearance of only a single peak around 24 h 
was predicted, as opposed to the double peak structure 
that we have observed. 

Spectral leakage can be illustrated with the follow- 
ing simple example, where the data contain a signal 
with a period of 20 d (frequency = 0.05 cycles d _1 ). 
If we sample this signal regularly once-per day (i.e. dt 
= I d), we will obtain the set of measurements indi- 
cated by the data points (dots) in Figure These 
measurements could also be produced with (at least) 
two other sine curves, with frequencies of (1 + 0.05) 
and (1 - 0.05) cycles d" 1 . So, if we were to calculate 
the periodogram we would get spurious peaks at 1/dt 
± 0.05 cycles d _1 . In the evenly sampled case one 
would only ever extend the power spectrum up to the 
Nyquist frequency of 0.5 cycles d _1 (i.e. 0.5x(f/dt) 
cycles d" 1 ), so this spectral leakage would never actu- 
ally be observed. 
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Figure 7: Three sine waves with frequencies of 0.05 
d" 1 (solid line), (1 + 0.05) d~ x (dotted line) and 
(1 — 0.05) d _1 (dashed line), sampled at 1 d time 
intervals. All three sine waves can be reproduced 
by the observations indicated by the solid dots, 
demonstrating the effect of spectral leakage 



The 'pseudowindow' can be used to determine the 
response of a given data analysis system to a single 
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frequency sine wave (Scargle 1982). All spectral leak- 
age effects present in a given data set will manifest 
directly in the pseudo-window, indicating where any 
enforced quasi-periodicities lie in the sampling rate. 
Pseudowindows can be generated by simulating a sin- 
gle frequency sine wave using the same sampling rate 
as the source of interest. The frequency of the sine 
wave should be high enough so that spurious peaks 
from negative frequencies do not overlap with those 
from positive frequencies. A frequency of 2xl0~ 4 Hz 
should be sufficient for most RXTE ASM data sets. 
The resulting Lomb-Scargle power spectrum of this 
simulated data set is the pseudowindow, and indicates 
where power will shift to due to spectral leakage. By 
translating the simulated power spectrum along the 
frequency axis so that the peak with the highest power 
(i.e. 2xl0" 4 Hz) lies at zero frequency, we can di- 
rectly determine the effective sampling rates that are 
causing the observed spectral leakage. This observed 
power spectrum is thus the convolution between the 
true power spectrum and the pseudowindow. Figure [H] 
shows pseudowindows generated from the data on GX 
339-4 and EXO 2030+375. Note that the power spec- 
tra are not completely symmetrical around zero fre- 
quency. After the zero frequency, the strongest power 
lies at 1.1604xl0~ 5 Hz (~0.997 d) for GX 339-4 and - 
1.18275X 10~ 5 Hz (~0.979 d) for EXO 2030+375, which 
are attributed to enforced ~24 hr periodicities in the 
light curves (see Table 1 for the a list of the enforced 
~24 hr periods for all six test sources). The ~96 
minute orbital period of RXTE is also manifest in both 
power spectra at ~±1.735xl0 -4 Hz. 



± the 'effective' sampling frequency (1/dt). 
spurious peaks (_F s i , -F s 2) will be seen at: 



F 



si 



= F in -(l/dt) 



^s2 = 



F in + (1/dt) 



So the 

(1) 
(2) 



For the case where the input frequency is less than 
the effective sampling frequency F g ^ is negative and so 
does not explain the appearance of two peaks within 
our frequency range, leading to the conclusion that we 
would only see F s 2- Taking into consideration that 
_F m also produces power at 
spurious peaks to consider: 



-F m , there are two more 



F 



s3 



= -i* - (1/dt) 



F 



s4 = 



-F m + (V*) 



(3) 
(4) 



As F g 3 will always be negative and thus not present in 
our frequency range, Equation (3) can be disregarded. 
That leaves three possible frequencies (F s ^, F s 2 or 
F34) that may be present. For the case where the 
input frequency is less than (1/dt), F s 2 and F s ^ will 
be present in the power spectrum. As _F m approaches 
(1/dt), the two spurious peaks will move apart. As 
F m approaches zero, the two peaks will move closer 
together, converging at (1/dt). For the case where the 
input frequency is greater than (1/dt), F s i and -F s 2 
will be present in the power spectrum. We compute 
the Lomb-Scargle periodogram for positive frequencies 
only so we can summarise Equations (1) to (4) in a sin- 
gle equation: 



Fs = |F in ±(lM)| 



(5) 



EXO 2030+375 



- km ■ 



0.0002 -0.00015 -0.0001 -0.00005 0.00005 0.0001 
FREQUENCY (HZ) 



00015 0.0002 



Figure 8: Pseudowindows generated from the data 
on GX 339-4 (lower panel) and EXO 2030+375 
(upper panel) . Power from the RXTE orbital pe- 
riod of 96 min is marked by a (positive) and a' 
(negative). The peaks due to the enforced ^24 hr 
effective sampling rate are marked by f3 (positive) 
and P' (negative) 



Figure [5] shows the periods at which we expect to see 
spurious peaks, as a function of the input period. 




1 10 
INPUT PERIOD (D) 



Figure 9: The periods at which we expect to see 
spurious peaks, as a function of the input period, 
obtained from Equation (5) with a sampling period 
of 1 d 



Pseudowindows are useful because they indicate 
where the spurious power will be seen for a certain 
sampling pattern. In the case of the RXTE ASM data, 
we will see spurious peaks at the input frequency (F m ) 



5.3 Results 

The effective sampling frequency 1/dt can be deter- 
mined directly from a pseudowindow generated from 
the data on each object. This value of 1/dt can then 
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be used in conjunction with Equation (5) to deter- 
mine the frequency of the spurious peaks at ~24 h in 
a given power spectrum with known period(s). Using 
this method, we can test whether the observed spu- 
rious 24 h periods are in fact a product of spectral 
leakage as described above. 

To test this hypothesis, we simulated a 2xl0 -4 Hz 
sine wave sampled at the same rate as each of our test 
sources, in order to determine the effective sampling 
rate 1/dt. Figure [TU1 shows pseudowindows between 
0.95 - 1.05 d for all six test sources, while the effective 
sampling periods are listed in Table Q with the uncer- 
tainties determined by taking the FWHM of the ~1 d 
peaks. To compare these sampling rates with the ac- 
tual spectral leakage peaks, we simulated 1 000 d sine 
waves with the same sampling rates. The resulting 
power spectra are shown in Figure ITTIfor between 0.95 
- 1.05 d. The effective sampling period peaks shown in 
Figure Hoi lie neatly between the two spectral leakage 
peaks for each source shown in Figure lTTI An interest- 
ing correlation between the location of the source and 
the sampling rate can be seen. Sources A, D, and F 
are all located in the northern hemisphere (i.e. pos- 
itive declination) and all have an effective sampling 
period of ~0.979 d. In contrast, sources B, C, and E 
are all located in the southern hemisphere (negative 
declination) and all have an effective sampling period 
of ~0.997 d. This apparent correlation is likely a result 
of slight differences in the semi-regularity imposed on 
northern and southern sources, possibly due to factors 
such as solar X-ray contamination, Earth occultation, 
and passage of the SAA. As the ASM is either turned 
off or the data discarded during each of these instances, 
the process inherently enforces semi-regularity on the 
sampling rates, possibly differing with source declina- 
tion. 



Table 1: Effective sampling 


rates dt & source 


dec 


Source 


dec 


dt 






(°) 


(d) 




Crab Pulsar 


+ 22.01 


0.9787 ± 3x10- 


-4 


LMC X-4 


- 66.37 


0.9973 ± 2x10- 


-4 


Vel X-l 


- 40.60 


0.9973 ± 3x10- 


-4 


Her X-l 


+ 35.34 


0.9785 ± 5x10- 


-4 


GX 339-4 


- 48.79 


0.9973 ± 2x10- 


-4 


EXO 2030+375 


+ 37.64 


0.9785 ± 6x10- 


-4 



The respective input periods were calculated for all 
observed ~24 h peaks with > 99% significance using 
Equations (2) & (4) and the effective sampling rates. 
All input periods that were present in the power spec- 
tra of Figure|5|were then simulated in light curves with 
the same sampling as the test sources. For all sources 
except the Crab pulsar and GX 339-4, these periods 
correspond either to the reported orbital and super- 
orbital periods of the sources, or to fractions of these 
periods. The periods used in the simulations for these 
sources are as follows: LMC X-4 - 30 d (super-orbital 
period); Vela X-l - 9.8 d (orbital period), 4.9 d (1/2 
orbital period); Her X-l - 35 d (super-orbital period), 
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Figure 10: Pseudowindows for the six test sources 
showing the ~1 d effective sampling rate dt. 
The apparent similarity of the pseudowindows for 
sources A/D/F and B/C/E appear to be a result 
of the sign of the declination of the sources (i.e. 
northern or southern hemisphere) 



17 d (1/2 super-orbital period), 1.7 d (orbital period); 
EXO 2030+375 - 46 d (orbital period), 23 d (1/2 or- 
bital period), 16 d (1/3 orbital period). As mentioned 
previously, the cause of the long-term periods in the 
Crab light curve is unknown, although the two periods 
used for the simulations (2644 d and 51 d) are possibly 
related to the solar cycle and precession of the RXTE 
orbit. The six periods used for the simulation of GX 
339-4 are due to two large scale flaring episodes in the 
light curve. These state transitions introduce a num- 
ber of trends and periods into the light curve and result 
in significant spectral leakage around 24 hr. 

The resulting power spectra are shown in Figure lT2l 
for periods between 0.8 - 1.2 d, with the real power 
spectra of Figure El shown for comparison. Spectral 
leakage can be seen to account for the most significant 
peaks around 24 h in each of the test source power 
spectra. Any periods above the 99% significance level 
that are not present in the simulated power spectra of 
Figure [T2l do not result from the leakage of power from 
long-term periods present in the real power spectra 
by the ~1 d effective sampling period. These periods 
are potentially real, although further analysis is re- 
quired before a firm conclusion can be drawn. Spectral 
leakage explains why the values of the ^24 h periods 
change between sources, and why the spurious period 
is not present in data from all sources. Each source 
monitored by the ASM has a unique sampling pattern 
so that while sources close to each other (as viewed 
from the Earth) have similar sampling patterns, no two 
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Figure 11: Power spectra of simulated 1 000 d 
sine waves with the same sampling as the six test 
sources. The double peaks are located either side 
of the effective sampling period of each light curve 
as shown in Figure [Till 



Figure 12: Simulated power spectra (grey lines) 
of the most significant periods present in the light 
curves of the six test sources with the real power 
spectra from Figure[3]shown for comparison (black 
lines). The dashed lines indicate the 99% signifi- 
cance levels for the real data 



are identical. While it is likely that ASM data for a 
given source have an enforced semi-regularity and thus 
an effective sampling rate of ~1 d, it is by no means 
guaranteed. In order for a significant spurious peak 
at ~24 h to be present in a light curve two conditions 
need to be met - the sampling rate must have a high 
enough degree of enforced regularity, and significant 
real low-frequency modulation(s) must be present. If 
one or the other of these conditions is not satisfied in 
a given light curve, the spurious 24 h peaks will not 
appear. 

6 Conclusions 

In depth analysis of the spurious 24 h periods seen 
in RXTE ASM data sets leads us to the conclusion 
that the observed periods are not dominantly a re- 
sult of daily variations in background levels. While 
events such as activation due to SAA passage may be 
a contributing factor, this is not the primary mech- 
anism. We instead conclude that spectral leakage of 
low-frequency modulation due to an effective sampling 
rate of ~1 d is the dominant cause of the spurious sig- 
nals. 

In order to identify a suspect period around 24 h as 
spurious, a simple technique can be employed. Accord- 
ing to Equation (5), re-binning (i.e. changing 1/dt) the 
light curve into sample sizes less than 1 d (e.g. 0.3 d 
averaged observations) should remove all peaks around 
24 h due to spectral leakage of any real periods greater 
than ~0.5 d, leaving behind only real periods (Dieters 



et al. 2005 employed this technique on GX 339-4 in 
an attempt to identify the reported ~1 d orbital pe- 
riod). To test this theory we rebinned the light curves 
of all six test sources using a bin size of 0.3 d. Power 
spectra of these rebinned light curves are shown in Fig- 
ure ED The significant peaks at ~24 h in each of the 
power spectra were shifted to 0.3 d, leaving no signifi- 
cant peaks behind. As a further check, we repeated our 
analysis using a bin size of 0.5 d and found that no sig- 
nificant peaks were present around ~24 h. In addition, 
any real periods around 24 h should produce spurious 
peaks of their own which can be identified using Equa- 
tion (5) and by determining the effective sampling rate 
through calculating the pseudowindow. These tech- 
niques will not distinguish a real period from one pos- 
sibly produced through daily variations in background 
levels. However, any spurious periods produced by this 
effect should have a peak very close to the rotational 
period of the Earth (i.e. 23.9345 h). Any suspect pe- 
riod not located at this exact value and not a product 
of spectral leakage could be considered to be real. Fur- 
thermore, filtering out all observations influenced by 
SAA activation should remove any power due to daily 
variations in background levels, leaving only the real 
signals and those due to spectral leakage. Additional 
evidence in the form of periodic variations at other 
wavelengths or data from other instruments/spacecraft 
missions would be advised before a firm conclusion is 
drawn. 
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Figure 13: Power spectra of the six test source 
light curves after re-binning using a bin size of 0.3 
d. No peaks above the 99% significance levels (in- 
dicated by the dashed lines) appear between 0.8 - 
1.2 d 
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